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A theoretical study of the magnon-drag Peltier and Seebeck effects in ferromagnetic metals is presented. A magnon 
heat current is described perturbatively from the microscopic viewpoint with respect to electron-magnon interactions 
and the electric field. Then, the magnon-drag Peltier coefficient IT mas , is obtained as the ratio between the magnon heat 
current and the electric charge current. We show that n mag = C mag T 12 at a low temperature T; that the coefficient 
Cmag ls proportional to the spin polarization P of the electric conductivity; and that P > for C raag < 0, but P < 
for C raag > 0. From experimental results for magnon-drag Peltier effects, we estimate that the strength of the electron- 
magnon interaction is about 0.3 eVA 3 ' 2 for permalloy. 

KEYWORDS: Peltier effect, Seebeck effect, magnon drag, spin caloritronics, spintronics, spin current, electron- 
magnon interaction 
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Understanding interactions between spin dynamics and 
transport phenomena in ferromagnetic materials is one of the 
central issues in spintronics and spin caloritronics. 1 ' In partic- 
ular, it is important to clarify the thermal properties of spin dy- 
namics, that is, magnon scattering mechanisms, at finite tem- 
perature. The magnon-drag thermoelectric effects proposed 
by Bailyn 2 ' are good subjects for investigating the thermal 
properties and transport phenomena of magnons because the 
effects are determined by substance-specific properties 3 ' and 
are sensitive to external magnetic fields. Moreover, the effects 
are useful in the study of fundamental electron-magnon in- 
teractions. Blatt et al. A) reported that the Seebeck coefficient 
S for iron takes a maximum value near 200 K, that S is un- 
responsive to doping with heavy atoms and annealing, and 
that S is well fitted by the form S = C\T + C 2 T 3 ' 2 (where 
C\ and C2 are constants). These differ from the features ex- 
pected for phonon-drag Seebeck effects 3 ' both quantitatively 
and qualitatively. Grannemann and Berger 5 ' phenomenolog- 
ically derived a convenient expression for the magnon-drag 
Peltier coefficient in ferromagnetic metals, assuming that the 
drift velocity of magnons, v mag , is proportional to the drift ve- 
locity of electrons (i.e., v mag = rjv e ), 
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where D is the spin-wave stiffness constant, e > is the el- 
ementary charge, n e is the number of conduction electrons 
per unit volume, V is the volume of the system, n(x) :- 
[expix/ksT) - is the Bose-Einstein distribution func- 
tion, Icb is the Boltzmann constant, T is the absolute temper- 
ature, and u> q is the energy dispersion relation for magnons. 
By using this expression, they found that at low temperature, 
Ilmag 00 T 5 ^ 2 , and that consequently, the magnon-drag See- 
beck coefficient S m ae = n'^g/T oc T 3 ^ 2 . In addition, by fitting 
eq. (1) to experimental data, they obtained the values 77 = 2.98 
for Ni69Fe3i and 77 = 2.20 for Ni66Cu34. Recently, Costache 
et al. 6 ^ have developed a magnon-drag thermopile that can 
cancel out all contributions except for that of magnons to 



thermopower, allowing the electric voltage induced by the 
magnon-drag Seebeck effect to be measured directly. They 
found that the Seebeck coefficient of a magnon-drag ther- 
mopile made of permalloy takes a peak near 180 K and is well 
fitted to IT^ag/T at temperatures below 180 K (77 = 3). As 
above, the phenomenological expression (1) can accurately 
describe the experimental results. However, there has been no 
progress on the theoretical front for understanding magnon- 
drag thermoelectric effects from the microscopic viewpoint. 

In this letter, we describe magnon-drag thermoelectric phe- 
nomena in ferromagnetic metals at the microscopic level and 
aim to obtain detailed information on the magnon-drag ther- 
moelectric coefficients. We calculate a magnon thermal cur- 
rent driven by an electric field using the Keldysh Green 
function technique. The mathematical form of the resultant 
magnon-drag Peltier coefficient conforms to Grannemann and 
Berger's expression (1) and a microscopic expression for 77 is 
obtained. Furthermore, we show that 77, a physically important 
feature, is proportional to the spin polarization of the electric 
conductivity. 

Let us consider a Hamiltonian describing ferromagnetic 
metals <H := % ka . s k<T c l ka _c k(T + £ 9 (o q b\b q +<H em , where c ka is 
an operator that annihilates the <x spin electron with an wave 
vector k, s kcr := h 2 k 2 /2m - <xA/2 -fi is the energy dispersion 
relation for electrons calculated from the chemical potential fi, 
m is the mass of an electron, A stands for the exchange split- 
ting, b q is an operator that annihilates a magnon with an wave 
vector q, and < K em is an electron-magnon interaction defined 
by 7 - 9 ' 
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b q C kfk+q{ + ^9 C k+qi C k\ 



where / represents the strength of the electron-magnon in- 
teraction. The interaction with a static homogeneous elec- 
tric field E := -dA(f)/df is given by 'H'(f) := -j e ■ 
A(t)V in terms of the electric charge current density operator 
j e := —(eti/mV)J^ kcr kclc,. We define the space-averaged 



*E-mail address: dmiura@solid.apph.tohoku.ac.jp 



1 



J. Phys. Soc. Jpn. 



LETTERS 



magnon heat current density operator by 
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7mag« := — 2j qu q b{{t)b q (t), 
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where operators are in the Heisenberg representation with re- 
spect to ( H + ( H'(t). Using the lesser function 10, n) of magnons 
D q (t,f) :- -(i/K){b' q (t')b q (t)y the statistical average of 
y'magW can be represented by 

(jmag(O) = ^ J] f) ' 

where (• • •) is a statistical average in H + when W(0 is 
regarded as the nonequilibrium part. We perform the lowest- 
order calculation for // m ag(/)) with respect to both 'H em and 
'H'(t) as shown in Fig. 1. In energy space, we have 



CO-), 



where AD q (co+, coJ) corresponds to the diagrams shown in 
Fig. 1 and is given by 



/}D q (co + , coJ) :- 
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A(Q) ■ d q (co+)A q (co+, co-)d q (co-) 



A q (co + , a)-) 



dE 



8kt(E + )gki(E-)g k+q i(E + co) 



+ gki(E+)gki(E-)g k -tf(E - co) 



where co ± := co + Q/2, E ± : = E ± Q/2, A(Q) := 
f d/ e' n ^ fi A(f), and gka-(E) and af 9 («) are the unperturbed 
Keldysh Green functions of the <x spin electron with the 
wave vector k and the magnon with the wave vector 
q, respectively. Taking the lesser component, we obtain 
d q (co + )A q (co+, co^d^co-)^ = d r q (co + )A I q (co+, co-)d q (coJ) + 
d q (co + )A q (co + , co-)d q (co-) + d q (co + )A q (co+, co-)d' q (co-) where 
d q (co) = -n(co)[d' q (co) - d l q (co)], and d q (co) and d q (co) are the 



unperturbed advanced and retarded Green's functions of the 
magnon, respectively. The part representing the electron-hole 
pair is, in the first order of Q., calculated as 
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A q (co+,co-) = A\(co+,co-)*, 



A q (co+, «_) = i3 
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J> f dE 



df(E) 
dE 



f(E + co)p k} (E) 2 g\ +q[ (E + co) 



+ {l-f(E- L0)) Pki (E?g\_^{E - co) 



where f(x) := [exp(x/k B T) +1] 1 is the Fermi-Dirac distri- 
bution function, g r kcr (E) := [E - s k(T + iS]' 1 = g\JET, and 
p k(T (E) :- -3 g\ (E) I Ji . Note that A q (co, co) vanishes iden- 
tically when performing k summation and E integral. The 



factor -df(E)/dE acts on the Green's functions of the elec- 
tron approximately like the Dirac's delta function 6(E) un- 
der the condition 5 » k B T. This approximation is valid at 
temperatures below about 100 K because 5 is on the order of 
10 -2 eV in metals. Considering only the contribution of the cr 
spin electrons to the Fermi surface \k\ = | Arpo-I and the long- 
wavelength region t^k-pa • q/m(co - A) <K 1, we obtain 

„ 2nm o> - cri 
A*(co + ,co-)^qn— T t . l — , (2) 
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where ov : = (7re 2 H 3 /3m 2 V)J] k k 2 p ka -(Q) 2 is the spin- 
dependent electric conductivity, 121 and we use the identi- 
ties / dE df(E)/dE = -1 and f dE f(E + co)df(E)/dE = 
/ d£[l - f(E - co)]df(E)/dE = (d/dco)con(co). Using eqs. (2) 
and (3), we have 

z v 2i 2 dp , r 

(jmag(O) - Jc^-y CO q j ( 



dw3 
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(co - A + i<5) 2 



n(co)d l „(co) 2 + co — - — ^-d T n (co)dl(co) 
q dco q q 



where j c : = (o"| + cr{)E is the electric charge current and 
P := ((Tf— cr^)/(o"f +<tj,) is the spin polarization of the electric 
conductivity. As a result, we obtain the magnon-drag Peltier 

coefficient Il mag := (jmagO)) 7 A/c (r = *,:y,z). This result 
indicates that the magnon-drag Peltier coefficient is propor- 
tional to the spin polarization. To evaluate the co integral, we 
express the self-energy of the magnon in terms of the Gilbert 
damping constant a as 13) 

1 



cf q (co) 



= d q (co)*. 



Considering a low-energy region co q < k B T in the q summa- 
tion and assuming a «c 1 and A <sc 6 /a, we obtain 

4I 2 DP6A 



n, 



mag 



3ae(A 2 + 6 2 ) 2 V 



dn(co q ) 

' diOa 



(4) 



A comparison between eqs. (4) and (1) affords the following 
microscopic expression for the phenomenological parameter 
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Q<(A 2 +5 2 ) 2 ' 

Grannemann and Berger 5 ' have experimentally obtained r\ 
2.98 forNi 69 Fe 3 i using n e = 3.1 x lO^A^ 3 in eq. (1). Putting 
these values into eq. (5), we can estimate the strength of the 
electron-magnon interaction to be / = 0.3 eV-A 3 ^ 2 by using 
typical values for ferromagnetic metals: A = 0.1 eV, 14) P = 
0.5, 14) 5 = 0.01 eV, 15) and a = 0.01. 15) 

Equation (4) could also be useful for determining the sign 
of P by a comparison with the experimental data. When co q = 
Dq 2 + A gap , eq. (4) can be written in the form 
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x I dx x 3/1 



A gap \ exp(x + A g . dp /k B T) 
k B T) [exp(x + A gRp /k B T) - 1] 



2 



J. Phys. Soc. Jpn. 



LETTERS 



_. C T 5 ' 2 

Since C mag can be determined experimentally by fitting a 
function with respect to T to the temperature dependence of 
the Peltier coefficient, we can evaluate the sign of P, that is, 
P > for C mag < and P < for C mag > 0. This is analo- 
gous to the situation where the type of charge carriers in semi- 
conductors can be determined by Peltier or Seebeck measure- 
ments. 

Finally, we refer to a microscopic description for Seebeck 
effects. The magnon-drag Seebeck effect is the counter phe- 
nomenon of the magnon-drag Peltier effect. Therefore, it is 
expected that the two effects will satisfy Onsager's reciprocity 
relation, that is, the magnon-drag Seebeck coefficient must be 



given by S r 



rimag/7^ The first work to microscopically 



describe a response to a temperature gradient was that of Lut- 
tinger. 16) Following that approach, we can describe Seebeck 
effects by introducing a Hamiltonian J Arh(r)4>\Jr, t), where 
h(r) is a total Hamiltonian density operator and <p^(r, t) is a 
pseudoscalar potential. In other words, we calculate the re- 
sponse of the nonequilibrium electric charge current to the 
Hamiltonian, and then replace V0 L (r, t) with (VT)/T. How- 
ever, here we consider the Hamiltonian -j mag ■ A^(t)V instead 
of J Arh{r)4>iir, t) from an analogy with the fact that the prod- 
uct of the charge density and the scalar potential corresponds 
to the product of the charge current density and the vector 
potential. Thus, we regard 



/mag 



2D v-i i 
A L (0V = 2_t <l<» q b\b q ■ A L (f) 



as a Hamiltonian describing a static homogeneous temper- 
ature gradient, and we assume that dAiXO/df corresponds 
to (VT)/T as E is given by -dA(f)/df. This Hamiltonian 
has some advantages: First, the diffusion ladder arising from 
short-range impurity scattering vanishes for a static homo- 
geneous temperature gradient. Second, it is easy to see that 
the Seebeck coefficient is associated with a charge current- 
energy current correlation function 17 ' because -/ mag ■ Al(/)V 
has the same form as the Hamiltonian — / e ■ A(t)V that we used 
to account for the static homogeneous electric field. 

Considering the lowest-order contribution shown in Fig. 2, 
we obtain the magnon-drag electric charge current j mag as 



.^)Vf fdEAUQ) 
27rmV/4 J J 



gkt(E+)gk+qi(E + uS)d q (u)+)d q (a)_)g k] (E_) 



+ gki(E+)g k - q i(E - u>)d q {b)+)d q (u-)g k i{E-) 

Through similar calculations as above, we obtain 

dA L (?) 4I 2 DP5A ^ 2 d "(<V 



J mag 



dt 



~3ae(A 2 +6 2 ) 2 V^ q ^ c\u q 



As a consequence, we can confirm that Onsager's reciprocity 
relation S mdg = n mag /r is reproduced by using the relation 
dAiXO/df — > Vr/T and the definition of the Seebeck coeffi- 



cient j' mag =: -o-S mag V7\ 

In summary, we have shown that n mag 



C mag r 5/2 and 



n m a g /r = C mag r 3/2 at a low temperature T from the 
microscopic viewpoint, and that the coefficient C mag is pro- 
portional to the spin polarization P of the electric conductiv- 
ity. Moreover, from a comparison with experimental results 
for magnon-drag Peltier effects, we estimate that the strength 
of the electron-magnon interaction is about 0.3 eVA 3 '' 2 for 
permalloy. 
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Fig. 1. Diagrammatic representations of magnon-drag Peltier effect. The 
solid and wavy lines represent the Keldysh Green functions of the electron 
and magnon, respectively; the dashed line stands for the vector potential; and 
the cross stands for the magnon heat current vertex. 
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Fig. 2. Diagrammatic representations of magnon-drag Seebeck effect. The 
solid and wavy lines represent the Keldysh Green functions of the electron 
and magnon, respectively; the dashed line stands for the pseudovector poten- 
tial; and the cross stands for the electric charge current vertex. 
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